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1. Introduction 


In classical mechanics, the Hamilton-Jacobi formalism is used to arrive at the 
frequencies of periodic systems, through the action variable, without having to ob¬ 
tain the complete solution of the equation of motion [1], Analogously, quantum 
Hamilton-Jacobi (QHJ) formalism has been developed [2] where the quantum ac¬ 
tion variable yields the exact energy spectra of the bound state problems, avoiding 
an explicit solution of the corresponding Schrodinger equation. The QHJ equation 
is equivalent to the Schrodinger equation and is written in terms of the quantum 
momentum function (QMF), p{x, E), which is the logarithmic derivative of the wave 
function. The quantization of the energy E arises from a suitable contour integral 
of the QMF, equated to integral multiple of h. In this paper we apply this formal¬ 
ism as developed by Leacock and Padgett to obtain exact energy eigenvalues for 
potentials which exhibit supersymmetry (SUSY) and shape invariance (SI) [3]. It 
is well known that for potentials having the above two properties, one can End out 
the eigenvalues algebraically. 

Our motivation in carrying out this explicit computation has been (a) to gain a 
better understanding of the working of the QHJ formalism and (b) to see what kind 
of singularities are possible for the QMF which play a signihcant role, as will be seen 
later, in determining the eigenvalues. It is hoped that the singularity structure of 
QMF will enable us to analyse exactness or nonexactness of the well known SUSY 
WKB formula [3,4]. 

The quantum Hamilton-Jacobi equation for a given potential V(x) is given by 

f,dHV(x,E) , fdWix.E)]" „ > 

7 ^ 

= pI{x,E), ( 1 ) 

where we have set 2m = 1 and where, W{x, E) is the quantum Hamilton’s charac- 
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teristic function. The quantum momentum function p{x, E) is defined as 


p{x,E) = 


dW{x,E) 

dx 


The QMF is required to satisfy the boundary condition 


( 2 ) 


lim p{x, E) —> E). 


Here Pc{x, E) is the classical momentum function, and is dehned in such a way that 
it has the value +\JE — V{x), just below the branch cut which joins the two turning 
points. 

Explicitly, the QMF is related to the wavefunction by 


., 1 d'lp 

P = • 

tp ox 


For the n-th energy level, if) has n zeros between the classical turning points. These 
zeros of the wave function correspond to simple poles of the momentum function 
with residue given by h/i for each of these zeros. The location of these poles is 
energy dependent and we shall refer to them as the moving poles. Other poles 
whose location does not depend on energy will be refered to as the hxed poles of 
QMF. Let C* be a contour enclosing the moving poles between the classical turning 
points. The integral 


J{E) = ^ £p{x,E)dx, (3) 

called the quantum action variable is obviously J{E) is equal to nh. Thus 


J{E) = nh, 


( 4 ) 


when inverted for E, gives the exact energy eigenvalues. 

In this paper we show how the QHJ method gives exact eigenvalues for SUSY 
potentials. SUSY quantum mechanics (QM) has been studied extensively in the 
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last decade. For a review of SUSY in QM, the reader is referred to recent review 
article by Cooper, Khare and Sukhatme [3]. In SUSY QM one considers a pair of 
potentials defined by (2m = 1). 

V± = ± h—u}{x), (5) 

(JjJb 

where u{x) is called the super-potential. The corresponding Hamiltonians 

H±=p^ + V^{x), 


can be factorized as 


//_ = A^A ( 6 ) 

= AAf, ( 7 ) 

where 

A = h-^+uj{x) ( 8 ) 

ax 

= -h^+uj{x). (9) 

ax 

Whenever the function exp {—uj{y)dy) is square integrable, it represents the 
ground state wave function of if_. Let 'ipo denote this ground state wave function. 
Then we have 

The eigenvalues and eigenfunctions for the two Hamiltonians H± are related. Using 
the intertwining relations AH_ = H^A and A^H^ = H_A\ it can be easily shown 
that (apart from the ground state wavefunction satisfying Aip^'^ = 0). 

Hence and are the energy eigenvalues for iL_ and respectively. It 

has been known that in addition to SUSY, if a potential has a discrete reparam¬ 
eterization invariance called shape invarience, the corresponding eigenvalues and 
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eigenfunctions can be explicitly obtained algebraically. In the next section we give 
the details of the steps in QHJ formalism required to solve one of these potentials. 
In Sec. 3. we shall take up the solution of a class of SUSY potentials. 


2. Eckart Potential 

In case of SI potentials both U+ and U_ are of the same functional form, albeit with 
different values of the parameter. In the following, we will consider only the U_ and 
apply the QHJ formalism to obtain the corresponding eigenvalues. In this section 
we shall obtain the energy eigenvalues for bound states of the Eckart potential given 
by 

V (x) = + — + AiA — ah)coseQb?ax — 2Hcoth ax. (11) 

A‘^ 

{x > 0) 


The corresponding superpotential is given by 

u(x) = —Acothax + —, (12) 

where A and B are constants and the quantum Hamilton- Jacobi equation is given 
by 

p^(x, E) — ^ = E — A^ -— — AiA — aKicoseob?ax + 2Hcoth ax . (13) 

ox A^ 

A simple way to obtain the eigenvalues of Eckart hamiltonian is to make use of the 
transformation 

y = cothax. (14) 

Under this mapping p{x, E) becomes a function of y when x is expressed in terms 
of y. We shall continue to use p to denote the function of y so obtained. Thus 
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p{y,E) means the QMF p{x,E), expressed in terms of the variable y. This will be 
understood for all ehange of variables to be considerd below. 


The quantization condition then becomes 



( 15 ) 


where Ci is the image of the contour C of enclosing the turning point, under 
the mapping x ^ y = coth ax. The QHJ equation written in terms of y is given as 



Note that the above mapping y = coth ar introduces additional singularities in the 
integrand at yEl. The contour integral Jci is calculated by deforming the contour, 
(see hg 1.) so as to enclose all the singular points of the integrand. The poles 
of QMF can be located using available results [5] on complex zeroes of solutions of 
linear differential equations. They are easily found by inspection for all the problems 
of interest in this paper. For evaluation of the integral ([T^ consider the contour 
integral Jr^ for a circle Fi? of radius R which is taken to be large enough so that 
outside F/j, p{y) has no singularities. The singular points of the integrand are the 
poles at y = ±1 and the moving poles enclosed inside Ci. Therefore, we have 



1 f pdy 1 


/ 

Jr 


2na Itr 1 — ‘J.TTa 


Denoting the integrals in (^) as Jr^, /cu Ai A 2 (^8-1)) rewrite ([T7|) as 


( 18 ) 


Ivr Ici T -^71 T -^72- 


The contour integral Rr is calculated by one more change of variable to z = 1/y. 
In terms of the variable the integral becomes 



( 19 ) 
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pdz 


70 


= L 


70 


(20) 

( 21 ) 


where 70 is a small circle in the z-plane enclosing only one singular point z = 0. It 
is worth reminding that both the contours are in the anticlockwise direction and the 
singularity at j/ —> cxd is mapped to the singularity ai z = 0. 

Therefore (|T^, (|T^ and (|^) give 


nah = — 


pdy 


1 


pdz 
1 — z'^ 


pdy 

1 - 2/2 


pdy 

i-y\ 


( 22 ) 


27r Jci 1 2/2 z,/i -L •^ 71 - 1 - y 

The calculation of various integrals requires behavior of the momentum function 
near the singular points. In particular, we need the value of the function at 2 / = ±1 
and its residue at z = 0. These are calculated by substituting appropriate Taylor 
or Laurent series expansion of p in the QHJ equation and solving for the hrst few 
coefficients in the series expansion. For example, for calculation of the contour 
integral around y = ±1, the series expansion of p{y,E) around 2 / = ±1 is used. 
The integrand suggests that we need to calculate only the coefficient of the constant 
term. We illustrate this for computing the integrals for this purpose we at hrst 
expand p{y, E) as 

p( 2 /, E) = aQ + aiiy - 1) + 02 ( 1 / - 1)^ +. 


and substite the expansion in (IT). Comparing the terms independent of y on both 
sides gives 


ai = E-A^ 


52 

I2 


+ 2B 


(23) 


Let /do be the constant term in the expansion of p{y, E) in powers of ( 2 / + 1). Then 
/do is similarly determined and is given by 




(24) 
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Note that, there are two roots for oq and /3o, corresponding to the two signs of the 
sqnare of the right hand side in (p3Dand Similarly for obtaining the integral 

we need to compnte the residne of p{z, E) at z = 0. We expand p{z, E) as 


p{z, E) =hi/z + ao + aiz H- 


(25) 


and snbstitnte ([^) the QHJ eqnation written in terms of the variable z 


p 


iah 


dp 


= E-A^-^-A{A-ah) 


z'^ J dz 

Comparing the coefficients of Xjz^ on the two sides gives, 

—iahEi{ah — 2A) 

A = -- 


+ 


2B 


( 26 ) 


( 27 ) 


We select the correct root for ao,/3o and bi from 0. @ and (p^) by comparing 
them with the answer for E = 0. It is straight forward to obtain the valne of bi, 
ao and Pq for E = 0, hy recalling that for zero energy the QMF is related to the 
snperpotential by p{x, E = 0) = iuj{x). Writing the snperpotential uj{x) in terms of 
z, we get 

, ^ ^ B 

( 28 ) 

The residne of uj{z) at z = 0 is —A. Comparing this answer with valne of bi, we see 
that the correct choice of bi is given by bi = —iA for all E. Similarly one looks for 
the coefficient of expansion of w in powers of ?/ ± 1 and this is compared with the 
valnes of ckq and (3q respectively for i? = 0. It is fonnd that these roots have relative 
opposite signs and are given by 

«o = -^E-A^-^ + 2B (29) 

/?o = ^E-A^-^-2B (30) 

It may be noted that the correct sign of the residnes calculated can also be hxed 
by taking h —>■ 0 and looking at the behaviour of Pc near the point of interest. This 









procedure, as originally suggested by Leacock and Padgett, is a bit complicated. 
For SUSY potentials under consideration in this article we have found it useful to 
follow the alternate procedure as explained above. 

The contour integrals and are, therefore, computed to be 




and 


A 

a 


Thus the energy eigenvalues are obtained from 


Jci = nh 


A i(3o iao 

a 2a 2a 


( 31 ) 


which on further simplihcation gives 

+ U - - ("“ft + . (32) 

One can arrive at the above energy eigenvalues using other mappings. Use of a 
different mapping will be considered in the next section and we shall then show how 
to obtain the bound state eigenvalues for other SUSY potentials. 
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3. Other SUSY potentials 

In this section we will show how the QHJ method could be used for other SUSY 
potentials. We will use y = exp(iQ;a;) mapping for the SUSY potentials involving 
trigonometric functions. The remaining potentials involving hyperbolic functions 
y = exp (ax) will be used. The treatment of each of these cases runs parallel to the 
treatment given to the Eckart potential in the previous section except for a new point 
which requires special attention. The use of the mapping y = exp (ax) gives rise to 
extra energy dependent poles in p(x, E) in the non-classical region. Computation of 
the exact energy eigen-values requires the knowledge of an integral along a contour 
enclosing these energy dependent poles. In all the cases investigated in this paper 
this integral can be related to the integral around the contour which encloses poles 
in the classical region on the real axis. 

In this section at hrst we shall work out the eigenvalues for the Eckart potential 
again, using the mapping y = exp (ax), but concentrating only on the new points as 
compared to the treatment in the previous section. 

The corresponding super potential for the Eckart potential written in terms of the 
variable y = exp (ax) is 


‘^(v) 




(33) 


and the QHJ equation is 




4A(A - ah)y^ 2B{y'‘ + 1) 

(y^ - 1 )^ (y^ - 1 ) 


(34) 


Now the equation 


4A{A - ah)y^ 2B{y^ + 1) 

^ ( j ,2 - 1)2 + ( j ,2 - 1 ) 


0 . 


(35) 


has four solutions for the turning points. These are shown as A, B, A' B' (see Fig. 
2). The moving poles of p are (i) on the real axis and between A and B and (ii) 
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on the real axis between A' and B' . It should however be noted, that two of these 
turning points A' and B' are in the non-classical region. We note that the symmetry 
y —> —y in (^) interchanges A with A' and B with B' . This symmetry implies 

Ici = Ic2, ( 36 ) 

where Ci and C2 are contours enclosing A, B, and A', S' respectively as shown in 
Fig. 2. Next note that now ?/ = 0 is a pole of the integrand in the action integral 
(see (^) below). Also p{y) has poles at y = 1 and y = —1 because the right hand 
side of QHJ is singular at y = 1 and y = —1. Introducing Jr^ for the large circle Fr 
enclosing all the singular points we arrive at 

= Ici + Ic2 + -^7i + A 2 + As- (37) 

Here are integrals along contours 71,72,73 which enclose the singular 

points y = l^y = —1 and y = t) respectively. Therefore, using (|5BD, the quantization 
condition 

Ic, = / —p{y, E) = nh. ( 38 ) 

iTia Jci y 

becomes 

Itr -Y.K = 2 ^^- (39) 

p 

The rest of the calculation is same as in Sec. 2 . and one easily arrives at 

- („a A a)2 ■ 

The calculation of eigenvalues for other SUSY potentials proceeds in a similar fash¬ 
ion. The results are summarized in Table I and Table II. The range of the variable x 
is —00 < X < 00, unless indicated otherwise. Expressions listed in the coloumn four 
are are the values for a contour 7 ^ enclosing only the singular point indicated 
in the third column of the table. The value of for the pole at 00 stands for 
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the value Jr^. For the hxed poles at ?/ = 0 and oo the square roots in the residue 
are found to have relatively opposite signs, on comparision with the coefficients of 
corresponding terms in the expansion of uj{x) for E = 0. The eigen-values listed in 
the last column are obtained using (l3^) . 

4. Conclusion: 

In conclusion, we have explicitly worked out the eigenvalues for the SUSY potentials 
using the QHJ method. Apart from checking the correctness of the formalism, this 
exercise provides insight into solvability of these potentials. The main effort involved 
in use of this scheme lies in selecting the correct root for the residues needed. This 
problem was solved here by comparing the answers obtained from the QHJ for 
E = 0 with that obtained from the superpotential. In general the choice of the 
correct root for the residue can be made by using the boundary condition on QMF, 
viz., p ^ Pc in the limit h —0 and where the branch of Pc is selected in such a way 
that it corresponds to the positive sign just below the cut on the real axis joining 
the physical turning points. 

The Leacock-Padgett method is a powerful method for obtaining the eigenvalues 
analytically as well as numerically and can be applied to other potentials. It is 
worth pointing out that for the potentials considered in this article, the SUSY- 
WKB approximation gives exact answer. Since QHJ as shown here also gives exact 
answer, it is natural to enquire as to the relation between these two approaches.This 
will throw light on the question of exactness of the SUSY WKB formula. This 
investigation will be reported elsewhere. 
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Table 1 : Trignometric Potentials. 

For all the potentials listed in this table the mapping y = exp(iaa;) is used. 


Name of 

Potential 

Location of 


Eigen 

potential 


hxed poles 

value 



0 

Ve + a^ 


Scarf I 

-A^ + (A^ + + Aah)x 

i 

-A + B 

{A + nah)"^ 

(Trignometric) 

sec^ ax — B {2A + ah) x 



-A^ 


sec ax tan ax 

—i 

-{A + B) 



(—7r/2 < ax < 12) 

oo 

-VE + A^ 





0 

-^E + A^ - B^/A^ + 2iB 


Rosen- 

A{A — Q;h)cosec^ ax 

1 

A 

+ByA^ 

Morse-I 

-A^ + B^A^ 



A — [A + nah) 

(trignometric) 

+2B cot ax 

-1 

A 

—B‘^/\a + nah) 


(0 < ax < tt) 

oo 

Je + A^- ByA‘^ - 2iB 








Table 2 : Hyperbolic Potentials. 


For all the potentials listed in this table the mapping y = exp (ax) is used. 


Name of 
potential 

Potential 

Location of 
hxed poles 


Eigen 

value 



0 

-iVE - H2 


Scarf II 

+ {B^ - - Aan)x 

i 

W-A 


(hyperbolic) 

sech^ ax + B {2A + ah) x 



A^ — {A — nah)"^ 


sech ax tanh ax 

—i 

-W -A 




oo 

WE - H2 




0 

-i^E - A^ - B^/A^ + 2B 


Rosen - 

A^ + B^/A^ 

i 

-A 

H2 + B‘^/A^ 

Morse II 

—A(A + ah) sech^ ax 



— {A — nah)"^ 

(Hyperbolic) 

-h2Btanh ax 

—i 

-A 

-B'^IW-nahf 



oo 

i^E - H2 - - 2B 




0 

i^E - H2 - R2/H2 - 2B 


Eckart 

A^ + B^/A^ 

1 

A 

H2 + ^2/^2 

(Hyperbolic) 

+A{A — ah) cosech^ ax 



— {A + nah)'^ 


—2Rcoth ax 

-1 

A 

-By{A + nah)^ 


(x > 0) 






oo 

-We - H2 - B'^/A^ + 2B 





0 

-WE - H2 


Generalised 

H2 + {B^ + AW Aah)x 

1 

-A + B 


Poschl- 

cosech^ ax — B{2A + ah) x 



H2 — (H — nah)"^ 

Teller 

coth ax cosech ax 

-1 

-{A + B) 



(x > 0) 

oo 

WE - H2 
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